Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear

Numerical Methods for Multilattices

Among the efficient numerical methods based on atomistic models, the quasicontinuum (QC) method has attracted growing interest in recent years. The QC method was first developed for crystalline materials with Bravais lattice and was later extended to multilattices (Tadmor et al, 1999). Another existing numerical approach to modeling multilattices is homogenization. In the present paper we review the existing numerical methods for multilattices and propose another concurrent macro-to-micro method in the numerical homogenization framework. We give a unified mathematical formulation of the new and the existing methods and show their equivalence. We then consider extensions of the proposed method to time-dependent problems and to random materials.Comment: 31 page

Study of space environment fabrication and repair techniques Final technical summary report

Joining techniques for fabrication and repair of structures in extraterrestrial environmen

Field theory simulation of Abelian-Higgs cosmic string cusps

We have performed a lattice field theory simulation of cusps in Abelian-Higgs cosmic strings. The results are in accord with the theory that the portion of the strings which overlaps near the cusp is released as radiation. The radius of the string cores which must touch to produce the evaporation is approximately $r = 1$ in natural units. In general, the modifications to the string shape due to the cusp may produce many cusps later in the evolution of a string loop, but these later cusps will be much smaller in magnitude and more closely resemble kinks.Comment: 9 pages, RevTeX, 13 figures with eps

Absorbing boundary conditions for the Westervelt equation

The focus of this work is on the construction of a family of nonlinear absorbing boundary conditions for the Westervelt equation in one and two space dimensions. The principal ingredient used in the design of such conditions is pseudo-differential calculus. This approach enables to develop high order boundary conditions in a consistent way which are typically more accurate than their low order analogs. Under the hypothesis of small initial data, we establish local well-posedness for the Westervelt equation with the absorbing boundary conditions. The performed numerical experiments illustrate the efficiency of the proposed boundary conditions for different regimes of wave propagation

Matching Conditions in Atomistic-Continuum Modeling of Materials

A new class of matching condition between the atomistic and continuum regions is presented for the multi-scale modeling of crystals. They ensure the accurate passage of large scale information between the atomistic and continuum regions and at the same time minimize the reflection of phonons at the interface. These matching conditions can be made adaptive if we choose appropriate weight functions. Applications to dislocation dynamics and friction between two-dimensional atomically flat crystal surfaces are described.Comment: 6 pages, 4 figure

Numerical Modeling of Transient Wave Propagation for High Frequency NDT

Electromagnetic nondestructive testing (NDT) methods use frequencies ranging from low (dc) to high (microwave) frequencies [1]. Applications of numerical methods to model two- and three-dimensional low-frequency (dc or eddy current) nondestructive testing phenomena, where displacement currents can be omitted, were extensively examined, [2,3]. These are all interior boundary value problems. Finite element study of ultrasonic wave propagation and scattering in metals, which is again an interior boundary value problem, was recently reported in [4]. However, modeling of wave propagation for high-frequency NDT problems have not yet been attempted. Recently, finite difference methods in time domain have been successfully applied to solve transient electromagnetic wave propagation problems over the atmosphere and the ground [5], and time-dependent eddy current problems [6]. The method used here is a generalization of this work and is designed for numerical modeling of high-frequency electromagnetic wave propagation arising from nondestructive testing applications. The physical situation includes examination of the scattering effects by cracks inside a piece of material (especially dielectrics) or due to surface variations when the material is illuminated by a TM plane wave. This leads to an interface type problem

Fast Hybrid Algorithms For High Frequency Scattering

This paper deals with numerical methods for high frequency wave scattering. It introduces a new hybrid technique that couples a directional fast multipole method for a subsection of it scattering surface to an asymptotic formulation over the rest of the scattering domain. The directional fast multipole method is new and highly efficient for the solution of the boundary integral formulation of a general scattering problem but it requires at least a few unknowns per wavelength on the boundary. The asymptotic method that was introduced by Bruno and collaborators requires much fewer unknowns. Oil the other hand the scattered field must have a simple structure. Hybridization of these two methods retains their best properties for the solution of the full problem, Numerical examples are given for the solution of the Helmholtz equation in two space dimensions.Mathematic

On the well posedness of the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's field equations

We give a well posed initial value formulation of the Baumgarte-Shapiro-Shibata-Nakamura form of Einstein's equations with gauge conditions given by a Bona-Masso like slicing condition for the lapse and a frozen shift. This is achieved by introducing extra variables and recasting the evolution equations into a first order symmetric hyperbolic system. We also consider the presence of artificial boundaries and derive a set of boundary conditions that guarantee that the resulting initial-boundary value problem is well posed, though not necessarily compatible with the constraints. In the case of dynamical gauge conditions for the lapse and shift we obtain a class of evolution equations which are strongly hyperbolic and so yield well posed initial value formulations

A Novel Long Range Spin Chain and Planar N=4 Super Yang-Mills

We probe the long-range spin chain approach to planar N=4 gauge theory at high loop order. A recently employed hyperbolic spin chain invented by Inozemtsev is suitable for the SU(2) subsector of the state space up to three loops, but ceases to exhibit the conjectured thermodynamic scaling properties at higher orders. We indicate how this may be bypassed while nevertheless preserving integrability, and suggest the corresponding all-loop asymptotic Bethe ansatz. We also propose the local part of the all-loop gauge transfer matrix, leading to conjectures for the asymptotically exact formulae for all local commuting charges. The ansatz is finally shown to be related to a standard inhomogeneous spin chain. A comparison of our ansatz to semi-classical string theory uncovers a detailed, non-perturbative agreement between the corresponding expressions for the infinite tower of local charge densities. However, the respective Bethe equations differ slightly, and we end by refining and elaborating a previously proposed possible explanation for this disagreement.Comment: 48 pages, 1 figure. v2, further results added: discussion of the relationship to an inhomogeneous spin chain, normalization in sec 3 unified, v3: minor mistakes corrected, published versio